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We applied the Recurrent Variational Approach to the two-leg Hubbard ladder. At half-filling, our 
variational Ansatz was a generalization of the resonating valence bond state. At finite doping, hole 
pairs were allowed to move in the resonating valence bond background. The results obtained by the 
Recurrent Variational Approach were compared with results fronr Density Matrix Renormalization 
Group. 
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I. INTRODUCTION 

In the hope to get a better understanding of strongly 
interacting systems, there has been considerable interest 
in ladder systems.! These ladder systems have proven 
to be a theoretical wonderland, both analyticallya and 
numerically S However, much of the analytic work done 
on ladders has been in weak coupling (or perturbatively 
in some parameter), namely because there are very few 
analytic methods at strong coupling. Exact diagonaliza- 
tion, Monte Carlo, and Density Matrix Renormalization 
Group methods have been the primary tools for studying 
these systems at strong coupling. Each of these methods 
has both strengths and weaknesses when considering the 
lattice sizes, temperatures, and couplings on can con- 
sider. 

With the ability to fabricate these materials,!] lad- 
ders are not just a theoretical playground. For ex- 
ample, in {yO^P^O-;, thece are well separated two-leg 
ladders composed of yC^Jj Also the cuprate-like mate- 
rial SrCu-iOz consists of weakly coupled CuC>2 two- leg 
laddersjj and the material SV2QU3O5 consists of weakly 
coupled CuO-i three-leg laddersE 

Recently, a powerful analytic method was developed 
to deal with strongly coupled quasi-one dimensional svfr. 
terns — the Recurrent Variational Approach (RVA).BEI 
This method is similar .in spirit to Wilson's Numerical 
Renormalization Grouplj and Wiiite's Density Matrix 
Renormalization Group (DMRG).liil The key idea in all 
of these methods is to build up the system by adding on 
sites at the boundary. However, the real power of the 
RVA is that, though analytic, the physics of the problem 
is taken into account in an unbiased way and elucidated 
quite clearly. For example, the importance of different 
configurations in the ground state wave function is deter- 
mined without any outside assumptions, and the physics 
of these configurations is made clear. 

In this work, we apply the RVA to the two-leg Hubbard 
ladder at strong coupling and small dopings. Though a 
considerable amount of work has been done on the Hub- 
bard ladder and many of its properties are known, we 
are unaware of any work which has put this information 



together and constructed a ground state wave function. 
Our goal in this work is to provide a simple physical pic- 
ture of what the ground state might look like and to go 
ahead and construct a ground state wave function. 

The Hamiltonian of the two-leg Hubbard ladder is 
given by 



H = -*E ( c L c i+*,* + h - c ) - ^ E ( c L c i+y, s + h - 

+U^2n iA n id , (1) 

i 

where c\ creates a fermion at site i with spin s, ni >s — 

c l s c i s> * ^ s ^ ne hc-ppmg matrix element along the chain, 
t± is the hopping matrix element perpendicular to the 
chain (i.e., along the rung), and U is the on-site Coulomb 
repulsion. Site i has coordinates (x, y) with 1 < x < N 
and y = 1, 2. It will also be convenient to introduce the 
following two operators: 



4; 



t J 






34 



c j,T c i4 ' 



which creates a singlet across sites i and j , and 
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c i,T U 



(2) 



(3) 



which creates a doubly occupied site. 

The rest of the paper is organized as follows. In Sec. II 
we consider the half-filled case. In Sec. Ill we consider the 
Hubbard ladder at small dopings. In Sec. IV we present 
our results for the ground state energies and compare 
with DMRG. Finally, in Sec. V we summarize our results 
and present some concluding remarks. 



II. HALF-FILLED HUBBARD LADDER 

We begin with the half-filled ladder. What are the in- 
gredients necessary to construct a wave function which 
captures the physics of the half-filled Hubbard ladder? A 
trail of clues has been laid by previous works. First of 
all, we know that at strong coupling the half-filled 
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FIG. 1. Ground state wave function for the 2x2 half- filled 
plaquette. A circled link represents a singlet bond. (See text 
for a full explanation of each state.) 



TABLE I. Values of the parameters for the half- filled 2x2 
plaquette (with t — tx — 1) which give the exact groundstate. 
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16 
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-a 3 


-a 5 
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-1.0 
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-0.1229 


-0,3 


-05 



Hubbard model is equivalent to the Heisenberg model. 
Secondly, it is known that the two-leg Heisenberg lad- 
der (and Hubbard ladder) have a spin gap and short 
range spin correlations. Furthermore, it has been shown 
that the resonating valence bond (RVB) state capti. 
the essential physics of the two-leg Heisenberg ladder £2 
Hence, the RVB picture should capture the essential 
physics of the half-filled Hubbard ladder at strong cou- 
pling. Thus we propose a generalization of the RVB state 
as our variational ansatz. 

A key property of the RVB state is that it is the exact 
solution to the 2x2 Heisenberg plaquette. Therefore, we 
base our generalized RVB state on the exact solution to 
the 2x2 Hubbard plaquette. Since the 2x2 plaquette 
will serve as the basis for our ansatz, we discuss it in 
detail below. 



A. The 2x2 Plaquette 



The 


ground state of the 2x2 plaquette is given by 
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I <P*) = A[l > l),(a i l) I> (l,2) I °> + A (l,l),(2,l)^(2,2) I °> 

+ A (l, 2), (2, 2)^(1,1) I °) + A (l, 2), (2, 2)^(2,1) I °) ' ( 5 ) 

I ipo) is shown schematically in Fig. YU. 

We will be mainly interested in the case t = t±. In 
Table H we list the values of the parameters for several 
values of U (with t = t± = 1). 



Notice that the solution to the 2x2 plaquette (with 
t = ij_ ) has D 4 symmetry. (D4 is the symmetry group of 
the square.) However, the ground state does not trans- 
form in the scalar representation of D4; it transforms in 
the B 2 representation of D4. (B 2 is the one-dimensional 
representation which changes sign upon rotation by 90° 
and reflection about the diagonals.) B 2 coincides with 
the standard d x 2_ y 2 symmetry. The B 2 representation 
is what forbids some configurations, as those shown in 
Fig. [2j, from appearing in the ground state wave func- 
tion. 

A few more words are in order about the ground state 
wave function. Notice that a\ = — 1 (i.e., the weight 
of the horizontal singlets is equal to (minus) the weight 
of the vertical singlets.) This is the RVB mechanism - 
the system lowers its energy by resonating between hori- 
zontal and vertical singlets; a\ is the "RVB parameter" . 
Also, notice that \ a 2 |, | 03 |= 0(af). This will play a 
role in constructing the Ansatz for the half- filled ladder. 



B. The Ladder 



Using the configurations of the 2x2 plaquette as the 
basis of our Ansatz for the ladder, a typical configuration 
for the ladder is shown in Fig. |3j. The ground state will 
be a superposition of all possible configurations of the 
type shown in Fig. |3j. Therefore, it seems like working 
with this kind of state will be a formidable task. Fortu- 
nately, the RVA gives us a straightforward way of deal- 
ing with such a state - generate it recursively. Specifi- 
cally, the RVA builds the ground state of a ladder with 
N + v rungs using the knowledge of the ground states 
of a ladder with N,N +1, ...,N + u-l rungs. This is 
achieved by recursion relations which express the ground 
state \ N + is) in terms of the ground states {| N + i)} 
with i = 0, . . . , v — l£l 

Using these ingredients, we consider the following 
ansatz for the half-filled Hubbard ladder which is shown 




FIG. 2. States which the B2 representation prevents from 
appearing in the ground state wave function of the half-filled 
plaquette. 



FIG. 3. Typical configuration appearing in the ground 
state wave function of the half-filled ladder. 




FIG. 4. The RVA ansatz for the half-filled Hubbard lad- 
der. Note that links or sites connected with a circle represent 
a singlet bond. 



N + 2) =| <j> ) N+2 | N + 1) + a | <f>i) N+2 \N + 1) 

+ P I </>2)n+1,N+2 I N) +7 I (/>3)N+l,N+2 I AT) 
+ £ I <MiV+l,iSr+2 I #) + V I 05)jV+l,iV+2 I AT) 

+ <5 I ^e)iv,JV+i,iv+2 I N - 1) 

+ £ | h)N,N+l,N+2 \N-1) (6) 



schematically in Fig. |J.t-3 For completeness, the states 
in terms of the operators in Eqs. @ and are given in 
Appendix A. 

A few words are in order about the configurations in 
our Ansatz. i) ft, the weight of the horizontal bond (rel- 
ative to the vertical bond) , is the "RVB parameter" . For 
a true RVB state, we would have ft = — 1. Although 
we do not expect ft = — 1, from work on the Heisenberg 
model the variational parameters were shown to evolve 
smoothly with system size;@ we expect | ft \— 0(1). ii) 
Suppose we iterate the recursion relations once. Then 
| <po) and | </>i) (and also | cfe)) generate the terms 
shown in Fig. Q In the ground state of the 2x2 pla- 
quette, the states in Fig. g(b) have weight 0(a), the 
states in Fig. 0(c) have weight 0(a 2 ), and the states in 
Fig. ||(d) do not appear. Since we expect the parameters 
to evolve smoothlylj wc expect i] sa —a 2 , (rj is behav- 
ing as a "counter-term"; it's job is to subtract off the a 2 
contribution from | <p\).) iii) Even though we no longer 
have D± symmetry, initial calculations showed that the 
states in Fig. ||(b) do not appear in the ground state 
of the Hubbard ladder. (This is another indication that 
the parameters evolve smoothly from the 2x2 case to the 
ladder.) Therefore, we do not consider them in what fol- 
lows, iv) Since the configurations | (j>e) n,n+i,n+2 an d 
I 4>t)n n+i,n+2 appear as intermediate states for the 
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FIG. 5. States generated by | 4>o) and | ^i) (and also | </>s)) 
by running the recursion relations once. 
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FIG. 6. Other "resonances" which play an important role 
in the RVB picture. 

resonances shown in Fig. 0, it is necessary to include these 
states in the Ansatz to give us an RVB state. This can 
easily be seen by comparing the weights of the coefficients 
as shown in Table || (see below). 

In order to compute the values of the coefficients, we 
treat them as variational parameters and minimize the 
ground state energy with respect to these parameters. 
The ground state energy and other quantities appear as 
recursion relations. It will be useful to define 

E N = (N\H N \N), 



D N 


- (N-l\ N {4> \H N \N), 


Cjsr 


- (TV - 1 U (0i \H N \N), 


Z N 


= (N\N), 


Y N 


= (N-l\ N {<h\N), 


Xjy 


= (N-1\ N (</>! | N) . 



They arc supplemented by the initial conditions 

Z = 1,Y = 0,X = 0, 
E = 0,D = 0,Co = 0. 



(7) 



(8) 



To determine the values for the variational parameters 
for a given (finite) value of N, we iterate the recursion re- 
lations and minimize the quantity En/Zn numerically. 
The actual recursion relations are quite unwieldly; we 
have relegated them, as well as their derivation, to Ap- 
pendix A. 

The values of the variational parameters for various 
values of U (with t = t± = 1) are shown in Table O. The 
results were obtained on a 2 x 32 ladder. Notice that i) 
| ft \— 0(1), and we have produced an RVB state, ii) 
| r] |= 0(a 2 ), and r\ is indeed behaving as a counterterm. 
iii) | S |, | e |w a/3. Therefore, these configurations are 
non-negligible, suggesting that the resonances shown in 
Fig. 11 are important to the RVB picture. 



III. THE DOPED HUBBARD LADDER 

Now, we consider the doped Hubbard ladder. In Ref. 9 
it was shown that hole pairs moving through an RVB 

TABLE II. Values of the variational parameters for a 2 x 32 
half-filled ladder with t — 1 1 = 1. 
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FIG. 7. Ground state wave function for the 2x2 plaquette 
with two holes. 



background captures the essential physics of the t — J 
ladder at small dopings. Since the t — J model is the 
large U limit of the Hubbard model, we expect this pic- 
ture to hold for the Hubbard model at large U and small 
dopings. Therefore, we consider an Ansatz of hole pairs 
moving through our generalized RVB background for the 
Hubbard ladder at small dopings. 

Since the structure of the hole pairs is based on the 
exact solution to the 2x2 plaquette with 2 holes, we con- 
sider the 2x2 case in detail below. 



TABLE III. Values of the parameters for the 2x2 plaque- 
tte with two holes (with t = t± = 1) which gives the exact 
ground state. 
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b 2 
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1.0 


1.2470 


0.3569 
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0.2100 


24 
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1.3420 


0.1483 



that the ground state of the half-filled ladder transforms 
in the B 2 representation, we see that the operator that 
creates a hole pair out of the undoped system has d x 2_ y 2 
symmctryO These facts continue to be true for low dop- 
ing in larger ladders. Also, notice that 63 = 0(a) where 
a is from the half- filled ladder. This will play a role in 
writing our Ansatz for the doped ladder. 



A. The 2x2 Plaquete with Two Holes 

The ground state of the 2x2 plaquette with 2 holes is 
given by 



\$) =\<Po)+h \^)+b 2 \^)+b 3 \ip$) 
where 
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I tpQ ) is shown schematically in Fig. m. 

In Table III we list the values of the parameters for 
several values of U (with t = t± = 1). 

Notice that the solution to the 2x2 case (with t = t±) 
has D4 symmetry. However, now the ground state trans- 
forms in the scalar representation of D4. Recalling 



B. The Ladder 

Using the configurations of the generalized RVB state 
as well as the hole pair configurations, a typical configu- 
ration for the doped ladder is shown in Fig. H. 




^ 



FIG. 8. Typical configuration appearing in the ground 
state wave function of the doped ladder. 

The ground state will be a superposition of all such con- 
figurations shown in Fig. @. Fortunately, we can generate 
such a state recursively. Specifically, we build the ground 
state of a ladder with N + v rungs and P + /1 holes us- 
ing the knowledge of the ground states of a ladder with 
N,N+l,...,N + v-l rungs and P, P + 1, . . . , P + fx 
holes. This is achieved by recursion relations which ex- 
press the ground state | N + v, P + fi) in terms of the 
ground states {| N + i, P + j)} with i = 0,...,v—l and 

j=0,...,/i.B 

Using the above ingredients, we consider the following 
Ansatz for the doped Hubbard ladder 



N + 2P+1) = I Mn+2 I N + 1 P + 1) + a \ <j> x ) N+ z \ N + 1 P + 1) + \ fa) N+i,Tr+2 I N P + 1) 
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FIG. 9. The RVA Ansatz for the doped Hubbard ladder. 

which is shown schematically in Fig. g. For completeness, 
the states in terms of the operators in Eqs. and H are 
given in Appendix B. 

A few words are in order about the configurations in 
our Ansatz. i) The Ansatz contains the configurations 
for the RVB state, as well as the configurations for holes 
bound in pairs. The hole pair states are based on the 
exact solution to the 2x2 plaquette with 2 holes, and 
our ansatz reproduces this exact solution, ii) The states 
I <l>a)N,N+i,N+2 and | 4>\)n,n+i,n+2 extend over three 
rungs. These states are necessary to allow the hole pairs 
to move smoothly through the RVB background, iii) 
Note that physically, the picture of hole pairs moving 
through an RVB background can only be appropriate for 
low dopings. 

To compute the values of the coefficients in our Ansatz, 
we treat them as variational parameters and minimize the 
ground-state energy with respect to these parameters. It 
will be useful to define 
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= (N-1P\ N . 

= (N-1P\ N . 

= (N P\N P) 

= (N-IP\ N . 

= (N-IP\ N < 

They are supplemented by the initial conditions 

Zn,p=n — 1, Yn,p=n — 0, A0v,p=at = 0, 
En,p=n = 0, Dn,p=n — 0, Cn,p=n = 0, 
F n <p,p = for F = Z, Y, X, E,D,C. 
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for given (finite) values of N and P, we iterate the recur- 
sion relations and we minimize the quantity En^p/Zn^p. 
The actual recursion relations are quite unwieldly; we 
have relegated them, as well as their derivation, to Ap- 
pendix B. 

What is the nature of the state we have constructed? 
In order to answer this question, we plot /?, A, and £ vs. 
doping. These parameters contain most of the physics of 
our Ansatz. f3 is the "RVB parameter"; A and C, are the 
weights of the hole pair configurations. The results were 
obtained on a 2 x 32 ladder. 

First, consider Fig. [Ty(a). (3 begins at 0{— 1) and in- 
creases (i.e., becomes less negative) with doping until a 
critical doping, x c , where it vanishes. Beyond this dop- 
ing, j3 is positive. This has also been found for the t — J 
ladder.B Upon doping, the hole pairs cause destructive 
interference which degrades the RVB mechanism. For 
x > x c , this destructive interference has driven (3 pos- 
itive, and it is no longer appropriate to think of our 
state as describing hole pairs moving through an RVB 
background E 

Similar to the t-J ladder, the difference between x < x c 
and x > x c can be attributed to two different internal 
structures of the hole pairs. For x < x c the hole pairs 
have a d x 2_ y 2 structure relative to the RVB background. 
For x > x c the hole pairs hawe an s-wave like symmetry 
relative to their background^ 

Now consider Figs. [K](b) and |Ty(c). First of all, notice 
that A > £. This shows the impoxtance of the diagonal 
frustrating bonds for all dopings.E£l Also, notice that A 
and C both reach their maximum at x = 1/2. At x = 1/2 
the system is essentially a large scale reproduction of the 
2x2 plaquette with 2 holesB Indeed, the values of A and 
C at x = 1/2 are similar to their values for the 2x2 pla- 
quette. 
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To determine the values for the variational parameters 
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FIG. 10. (a) 13 vs. doping, x, for a 2 x 32 ladder at U=8 
(solid line), U=12 (dashed line), and U=16 (dashed-dotted 
line, (b) Same as (a) for A. (c) Same as (a) for £. 



IV. GROUND STATE ENERGIES 



First, we show results for Energy per site vs. U in 
Fig. 111]. For comparison, DMRGj results are presented 
for the same set of parameters.!^ At half-filling, as we 
would expect, our Ansatz is most accurate for large U 
and large t±. The ground state energy per site for a 
2 x 32 half-filled ladder as a function of U for various t± 
is shown in Fig. |ll|(a). For U — 8 and t± = 1, the energy 
from the RVA agrees with DMRG to within 90% and 
improves as U or t± is increased. Up to about U = 10, 
longer bonds (extending over at least 3 rungs) are coming 
into play. These states should be included in the Ansatz 
to further improve the overlap with the ground state. At 
U — 16 and t±_ — 1, our ansatz gives a ground state 
energy within 94% of the DMRG result. 

It should be noted that for the Heisenberg ladder, the 
RVB state gives a ground state energy within 94% of true 
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FIG. 11. (a) Ground state energy per site at half-filling 
(x = 0) vs. U for £j_ = 1 (solid line), t± = 2 (dashed line), 
£_L = 3 (dashed-dotted line). For comparison, DMRG results 
(shown as x) are presented for the same set of parameters, 
(b) Same as (a) except for x — 1/8. (c) Same as (a) except 
for x — 1/4. 



ground state energy, obtained from DMRG.H A recent 
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FIG. 12. "Pair-breaking" configurations which are playing 
a rather large role in the ground state of the doped Hubbard 
ladder. 



DMRG study of different ladder models found that the 
Hubbard model and Heisenberg model begin to agree 
only for rather large U {U ~ 16). ED Therefore, it is not 
surprising that the RVB picture becomes as good for 
the Hubbard model as it is for the Heisenberg model at 
[7«16. 

Fig. |ll](b) shows the ground state energy of the RVA 
Ansatz as a funtion of U for various t±_ for a doping of 
x = 1/8 on a 2 x 32 ladder. Again, we show energies 
obtained from DMRG for the same set of parameters. 
For U — 16 and t± = 1, the two energies agree to only 
within 77% and improves slightly as t± is increased. For 
example, at U = 16 and tx = 2, the overlap of energies 
increases to 87%. Further discrepancies occur when the 
doping is increased to x = 1/4 (see Fig. |ll](c) ). 

The differences in ground state energies occur due to 
the importance of "pair-breaking" configurations, like 
those shown in Fig. H3. The weights of these types of 
states increase as we move away from half-filling. Conse- 
quently, to further improve the RVA Ansatz, such states 
must be included in the wave function. Note also that 
our Ansatz would be essentially exact for the case where 
hole pairs are well localized on a rung. For the t — J 
model with J rung ~^> J c haimt, pairs are well localized 
along a rung, and the ground state is essentially a prod- 
uct of rung singlets and rung hole pairs. However, for 
the Hubbard model at strong coupling (i.e., U 3> t±,t), 
this is not the case. Holes would always rather occupy 
adjacent rungs, even for t± » t, since this minimizes 
the Coulomb energy from doubly occupied sites. To see 
this consider a 2 x 2 plaquette with 2 holes; let t -C t± 
and t, £j_ <C U. With 1 particle on each rung (i.e., one 
hole on each rung), the ground state energy is approxi- 
mately — 2ij_; with both particles on the same rung (i.e., 
both holes on the same rung), the ground state energy 
is approximately — J = — At\/U. Therefore, at large U, 
the particles would rather occupy adjacent rungs. (See 
Fig. If) 

The situation we have with the doped ladder is similar 
to what we had for the half-filled ladder in the early 
stages of this work. We found that without the states 



(a) 



(b) 



V 



FIG. 13. For t << t±, hopping along the rung dominates. 
By putting the particles (or holes) on adjacent rungs, we can 
have the situation shown in (a). However, by putting both 
particles (or holes) on the same rung, we get the situation 
shown in (b), which is energetically unfavorable. 




FIG. 14. Ground state energy per site vs. doping for a 
2 x 32 ladder (with t — t x = l)for U - 8 (solid line) and 
U — 16 (dashed-dotted line). For comparison, DMRG results 
for U = 8 (x) and U — 16 (*) are shown. 

| 4>q) and | <f>-j) which extend over three rungs (see Fig. 0), 
the RVA did not accurately reproduce the ground state 
even at extremely large U . However, once we included 
| 4>e) and | fa), the results from the RVA improved dras- 
tically. Based on these results, we expect the RVA to 
greatly improve by including the states shown in Fig. |l2]. 
In Fig. O we plot Energy per site vs. Doping for a 
2 x 32 ladder for U = 8 and U = 16 (with t± = 1) in 
order to better understand the region of validity of our 
RVA ansatz. Again, we see good agreement with DMRG 
results at half-filling. However, as soon as we dope, con- 
figurations like those shown in Fig. [L2 are also important. 

It is interesting to note that the idea of hole pairs mov- 
ing through the RVB background seems to more accu- 
rately represent the ground state of the t — J model than 
the Hubbard model. Using the well known relation at 
strong coupling, J«|-, the RVA agrees to within 92% 
of the true ground state energy of the t — J ladder for 
J = 0.5 (U = 8) at a doping of x — 1/8. There are 
two ways to interpret this; either the t-J model supports 
pairing better than the Hubbard model, or we must view 
the hole pairs in the Hubbard model as having a larger 
size (i.e. larger coherence length.) 



V. CONCLUDING REMARKS 

To summarize, we applied the Recurrent Variational 
Approach to the two-leg Hubbard ladder. Our results 
were in qualitative agreement with previous results on 
the Heisenberg and t-J ladders. For the half-filled lad- 
der, the generalized RVB state became more accurate in 
the parameter regime where the Hubbard and Heisen- 
berg ladders were shown to coincide. However, compar- 
ison of the RVA with DMRG for the doped ladder indi- 
cate that hole-pairs moving through an RVB background 



is incomplete; it does not capture the essential physics. 
"Pair-breaking" configurations are also necessary to cap- 
ture the essential physics. 

As we saw, the strength of the RVA is the ease in which 
we could extract the physics. We were able to see the im- 
portance of the configurations in our Ansatz quite easily. 
Furthermore, the RVA has a natural way in which to in- 
clude longer bonds in the Ansatz to more accurately rep- 
resent the groundstate wavefunction. The importance of 
such additional states to the physics of the ladder is not 
easily probed with other techniques. 

Generalized RVB states similar to ours have been cop. 
sidered previously for the half- filled Hubbard ladder S3 
Fano et. al. were even able to produce an Ansatz coming 
within 98% of the true groundstate energy for a (half- 
filled) 2x4 ladder at U — 16. (Their Ansatz included 
diagonal bonds of length y/E.) However, none of these 
works considered the doped case. Using the approach 
in Ref. 17 (in terms of dimer coverings), it appears to 
be a formidable task to consider doping. This is one of 
the strengths of the RVA; doping is handled rather easily. 
Even though our results for the doped ladder showed that 
hole pairs moving through an RVB background is incom- 
plete, the RVA offers a straightforward way to improve 
the situation, namely include "pair-breaking" configura- 
tions (shown in Fig. 12 ) in the Ansatz. 



Another (and probably better) way to improve the sit- 
uation for tbe.|dpped ladder is to consider a Matrix Prod- 



uct Ansatz.1 



A Matrix Product Ansatz can be gen- 



erated by first order recursion relationsBO In the RVA, 
the size of the hole-pairs are fixed. (In our case, the 
hole pairs had a size of one lattice spacing.) However, 
by construction, the Matrix Product Ansatz takes into 
account hole pairs of arbitrary size. We leave this (and 
other possibilities) for future work. 
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APPENDIX A: THE HALF-FILLED LADDER 

The states in the Ansatz for the half-filled ladder of 
Eq. H are given by 



I 0O>AT+2 - A (7V+2,l),(7V+2,2) I °>W+2 , 

I <t>l)N+2 = D \ N+2 ,l) I 0)^+2 + D (N+2,2) I °)w+2 , 
(>2)N+l,N+2 = A (JV+l,l),(JV+2,l) A (JV+l,2),(jv+2,2) I 0)aT+1,JV+2 , 
h)N+l,N+2 = D (N+1,1) D (N+1,2) I ty N+l.N+2 + D (N +2,1) D {N+2,2) I ®) N+1,N 2 , 

^ r>t in\ xAt r>t 



I <Pi)n+1,N+2 = A (N+1,1).(N+2,1) D (N+1,2) I °) N+l.N+2 + A (N+1,1)XN+2.1) D (N+2.2) I °) N+l.N+2 
+ A (N+1,2)„(N+2,2) D (N+1,1) I °)n+1,N+2 + A (JV+1, 2), (JV+2, 2) D (N+2,1) I °)n+1.N+2 
I 4>h}n+1,N+2 = D (N+l,l) D (N+2,2) I ®)n+1,N+2 + D (N +1 ,2) D [N +2 ,1) I ®) N+l,N+2 , 



I 4>6}n,N + 1,N+2 — A (Af 1 



A 



D 



(N,l),(N+2,2) L - i (N+l,l),(N+2,l)^(N,2) 



0)n, 



N+l.N+2 



A 
A! 



A 



D 



(N,l),(N+2,2)'- i (N+l,l),(N+2,l)^(N+l,2) 



\0)n, 



N+l,N+2 



^A 



sD 



(N,2),{N+2,l) L ^{N+l,2),(N+2,2)^(NS) I ®)n,N+\,N+2 
+ A (JV,2),(A f +2,l) (JV+l,2),(JV+2 : 2)^(Ar+l,l) I 0)iV,JV+l,iV+2 : 

nt ] 0)n,i 



h)N,N+l,N+2 - A (N,1),(N+2,2) A (N,2)AN+1,2) D (N+1,1) 



,N+l,N+2 



A 
A 
A] 



A 



D 



(N,l),(N+2,2) L - i (N,2),(N+l,2)- L ^(N+2,l) 



A 



D 



(N,2),(N+2,l)'- i (N,l),(N+l,l)^(N+l,2) 



A 



,-D 



0)n, 
0)n, 

0)nj 



N+l.N+2 



N+l.N+2 



'(N,2),{N+2,l)^(N,l),{N+l,l)-^{N+2,2) I "/J*',N+l,N+2 ■ 

To derive the recursion relations, the following inner products are necessary: 



(1) 



N+2(4>0 I 0o)jV+2 = 2, at +2 (0i | 4>i)n+2 = 2, N+l,N+2(4>2 \ 4>2)n+1,N+2 = 4, 
N+l,N+2{4>3 | 4>3!N+l,N+2 = 2, N+1,N+2\4>A I 04)jV+l,iV+2 = 8, N+l.N+2^5 I 05)iV+l,iV+2 = 2, 

iV,iV+l,iV+2(<^6 | 06)iV,i\r+l,iV+2 = 16, N,N + l.N+2(4'7 \ &) N,N+l,N+2 = 16, 

(N + 1 \ N +2 (00 I 2 )iV+l,iV+2 | JV) = (-1)<JV + 1 | 0o)n+i | AT), 

(N + 1 | w+2 (0x I 5 )iv+i,Ar+2 I N) = (N + 1 I ^)^ +1 I AT), 

(N + 1 Iat+2 (00 I 0e)iV,iV+l,iV+2 I JV — 1> = (~l)(iV + 1 I (t>i)N,N+l \N-1), 

(N \n+i,n+2 (04 I <h)N,N+i,N+2 \N-1) = (-4)<iV I o )at I N - 1) . (2) 

Using these inner products, a straightforward calculation gives the following (coupled) recursion relations: 

Z N +2 - (2 + 2a 2 )Z N+1 - 2/3Y N+1 + 2ar,X N+1 + (4/3 2 + 2 7 2 + 8£ 2 + 2r?)Z N 

+ {16S 2 + 16e 2 - 16^)Ziv-i - 8&Y N + 88eY N _i , 

Fat+2 = 2Z N+1 - [3Y N+1 - 8$8Z N _ 1 + 4 < JeY 2V _ 1 , 

-X"iv+2 = 2aZ^v+i + r^Xjv+i , 

E N+2 = (2 + 2a 2 )^jv+i + (-8*j_a + 2[/a 2 )ZAr+i - 2(3D N+1 + 2ar)C N+1 

+ (8*£ + 4t ± af3)Y N+1 + (~4txV ~ 8ta£ + 2[/a?7)X i v + i 

+ (4/3 2 + 2 7 2 + 8£ 2 + 2r, 2 )E N + (-32*/3£ + 4[/ 7 2 - 16< 7 £ + 8U£, 2 - 16i£?/ + Wr, 2 )Z N 

+ (16t(38 - 8ie + 8ta 2 e + 16t/3e + 8i 7 e - 8U£s + 8ti]e - 8tS)Y N + 8taSX N 

+ {32t/3e + 16US 2 + 16Ue 2 + 32t/38 + 16£ 7 <5 - 16U£6 + 16tr]S + m±a6g)X N -i - 8&D N 

+ {16S 2 + 16e 2 - 168£)E N -i + 8fcL>Ar_i + (8USe - !U^a5e)Y N - 1 , 

D N+2 = 2E N+1 - 4t ± aZ N+1 - f3D N+1 + tt£Y N+l - 2t ±V X N+1 + (-Us - U8)Y N + tta8X N 

+ {16t(35 + 8tj5 - 8U^5 + 8tr]S + 16tps)Z N ^ - 8£8E N -i + 4feD„_i + 4U8eY n _ x , 

C N+2 = 2aE N+1 + (-4«_l + 2l7a)^ +1 + 2^/31^+1 + (-4i£ + C/??)X J v +1 + r/C w+ i 

+ 16ij.^iv_i - 8t ± 5£Y N -! + 4taeY N . (3) 



APPENDIX II: THE DOPED LADDER 

For the doped ladder (see Eq. [FI]), | 0o)jv +2 ,| 0i)jv+2, I 02)iV+i,iV+2 >| fa)N+i,N+2, I <Mjv+i,jv+2 ,| 0s);v+i,Ar+2, 
I 4 > g}n,n+i,n+2, and | 4 , y)n,n+i,n+2 are the same as the half-filled case, and 

1^+2 = 10)^+2, 

I 0l)jV+l,iV+2 = ^(JV+l,l),(JV+2,2) I 0)iV+l,JV+2 + A(AT+2,1),(JV+1,2) I ®) N+l.N+2 , 

I 02)jV+l,iV+2 = ^(JV+1,1),(JV+2,1) I 0)iV+l,iV+2 + ^(AT+l,2),(AT+2,2) I ty N+l,N+2 , 

I <f>3)N,N+l,N+2 = A (AT,l),(Ar+2,l) A (A'+l,2),(A'+2,2) I Q )n,N+1,N+2 + A (Ar,2),(Af+2,2) A (Af+l,l),(7V+2,l) I Q )n,N+1,N+2, 

I 04)iV,JV+l,JV+2 = A (AT,l),(Ar +2 ,l) A (A',2),(Ar+l,2) I °)n,N+1,N+2 + A (AT,2),(AT+2,2) A (Af,l),(A'+l,l) I °) N,N+l,N+2 ■ (1) 

To derive the recursion relations, we use the inner products from the half-filled case as well as the the following: 

Af+2(0o I 0o)jV+2 = 1) JV+1,JV+2(01 I 0l)jV+l,JV+2 = 4, N+l,N+2(4>2 I 4>2)n+1,N+2 = 4, 

N,N + l.N+2(4>3 I <p3/N,N+l,N+2 = 8, JV,JV+l,JV+2(04 I 04/JV,JV+l,iV+2 = 8, 

(N + 1P + 1 Iat+2 (0o I 03)^+1,^+2 | TV - IP) = {-l)(N + IP + 1 | <t>i)N,N+l \N-1P), 

(NP \ N +l,N+2 (01 I 4>i)N,N+l,N+2 | N - IP) = (-2)(7VP I o )at I N - IP) . (2) 

Using these inner products, a straightforward calculation gives the following (coupled) recursion relations: 

^v+ 2 ,p+i = (2 + 2a 2 )Z N+hP+1 - 2I3Y N+1 . P+1 + 2ar)X N+1>P+1 + (4/3 2 + 2 7 2 + 8^ 2 + 2t ? 2 )Z jv ,p+i 

+ (165 2 + 16e 2 - 16^)Zjv-i,p+i - 8^eFAr jP+1 + 8feY^_i,p+i + ^ 2 ^jv+i,p + (4A 2 + AC 2 )Z N . P 
- 4:\vY N:P + (8/i 2 + 8^ 2 - 8\n)Z N -i tP + 4nvY N -x,p , 

Fv+2,p+i = 2Zjy + i : p + i — /?Yjv + i,p+i — 8^<5Zat_i^p + i + 4<5eYjv-i,p+i — 4^A^Af_i ! p + 2/xu1jv-i,p , 



A"jV+2,P+l 
En+2,P+1 



2aZjv + i i p + i + i]Xn+i,p+i , 

(2 + 2a 2 )E N+ i tP+ i + (—8t±a + 2Ua 2 )Z N+ i^ P+ i - 2/3Djv+i,p+i + 2ar?Cjv+i,p+i 
+ (-U±T) - 8ta(, + 2Uarf)X N+ltP+ i 
2f] 2 )E N . P+1 + (-32i/3£ + 4U-f 2 - 16t 7 £ + 8U£, 2 - 16i£r? + Wr, 2 )Z N , P+1 
16t/3e + 8tje - 8U& + 8t-qe - 8tS)Y N , P+1 + 8ta8X NyP+1 
16US 2 + 16Ue 2 + m/36 + 16*7(5 - 16U& + 16tr]d + 32t±a6£)X N _ 1>P+1 - 8&D N , P+1 



(8t£ + 4t±ap)Y } 

(4/3 2 + 2 7 2 + 8<e 
{16t/3S - 8te 
(32i/3e 



8to 2 e - 



+ (16(5 2 + 16e 2 - 16^)^iv-i,p+i + 86eD N _ 1>P+1 + (8USe - 16t±a6e)Y N _ 1>P+1 

+ E N+1 , P + (4A 2 +4C 2 )^jv,p + (8m 2 + 8^ 2 -8A^)^jv-i,p - UXY N+1 , P - 4t(X N+hP 

+ (-16ij_AC - 8iA - 8taC,)Z N , P + {^v + 8t ± (v + 4i/x + 4H>)Y N ,p 

{-I6t(3v + 8tn + 16£_i_Cm + 16t±aA^ - 16t/3fi,)Z N -i tI 

Rt. i rvii.i/V*T i d . 



— 4\vDn, P + 4//i/Djv-i,p 



£> 



JV+2,P+1 — 



c 



iV+2,P+l 



- (-16i/3 
8ij_a^^yjv_i ,p , 

2-Ejv+i,p+i — 4ij_aZ]v+i,p+i — /3-Djv+i,p+i + 4££Yjv+i,p+i — 2£^?7AOv+i,p+i 
(-4te - 4M)F jv ,p+i + Ata8X NtP+l + (164/35 + 8i 7< 5 - 8E/£<5 + 8**7^ + 16t0e)Z N _ 1: 
8£6E N _ liP+1 + 4fcD„_ liP+ i + 4C/feF JV _i,p + i - UXZ NtP + (2tv + 2t/j,)Y N , P 
4A/i-Ejv-i,p + 2^zaDjv_i ! p + (42^ + 8t±(n)Z N _i jP , 

2a^iv + i,p+i + (-4t ± + 2L r Q !)Zjv+i,p+i + 2^/31^+1^+1 + (-4t£ + Urj)X N+ltP+1 
vCn+i,p+i + 16tj.^Zjv-i,p+i - 8tj_feFiv-i,p+i + 4tasY N ^ P+1 - 4t(Z NiP 
8tj_A/xZ/v-i.p — 4ij_/zz^yjv-i ! p . 



,p+i 



(3) 
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